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JL ■ Abstract 

Q ■ The elastic properties of rutile transition metal dioxides XO2 (X=Ru, Rh, Os, and Ir) are 

cj ! 

investigated using first-principles calculations based on density functional theory. Elastic constants, 
bulk modulus, shear modulus, and Young's modulus as well as Possion ratio are given. Os02 and 
ly-. Ir02 show strong incompressibility. The hardness estimated for these dioxides shows that they 

^f^ . are not superhard solids. The obtained Debye temperatures are comparative to those of transition 

L^ ! metal dinitrides or diborides. 
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^. • PACS numbers: 61.50.-f, 62.20.-x, 71.15.Mb, 71.20.-b, 77.74.Bw, 78.20.Ga 
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I. INTRODUCTION 

The rutile-type dioxides of transition metals (TMO2) possess a variety of electrical, me- 
chanical, magnetic, and optical properties. It's electronic structure and electrical proper- 
ties are widely studied due to their importance in technology. For example, OSO2 and 
Ru02 can be used as a possible electrode material in high-energy- density storage capaci- 
tors. Many of the TMO2S, such as RUO2, Ir02 and Ti02, present a rich phase diagram. 
They can transform to other structures under high-pressure conditions. The first theoreti- 
cal investigation on the electronic structures of the rutile transition metal dioxides (OSO2, 
Ir02, and RUO2), were reported by Mattheiss using a non-self-consistent angmented-plane- 
wave-linear-combinations-of-atomic orbitals method. [1|] Subsequently, plenty of theoretical 
researches, including on the electronic, optical and structural properties were performed for 
rutile TMO2. Xu et al studied the electronic and optical properties of Nb02, RUO2, and Ir02 
using self- consistent calculations based on the linear mufiin-tin-orbital method and atomic 
sphere app.oxi,na«„. fl KrasovsU . .1 calculated the optical a„d pKotoe.ect™ properties 
of Ru02 using the ah initio self-consistent energy-band structure. [3|] Recently, the electronic 
and optical properties of Ru02 and Ir02 were reported by first-principles self-consistent 
electronic structure calculations based on the full-potential linearized plane wave method. ^ 

In addition, the structural and mechanical properties of rutile transition metal dioxides 
have been widely investigated theoretically and experimentally. Structural phase transitions 
in Ir02|5[ and Ru02|6|] were reported based on first-principles calculations. Elastic proper- 
ties of several potential superhard Ru02 phases were studied by first-principlesplane-wave 
pseudopotential and full-potential linearized augmented plane-wave methods. [71] Yen et al 
measured the first-order raman specta of Os02 at room teniperature. 8| Recently, the lat- 
tice dynamics of RUO2 was given by theory and experiment. J9| However, to the best of our 
knowledge, few research on the elastic properties of Ir02, OSO2 and Rh02 has been reported 
so far. 

In this paper, we focus on investigating the elastic properties of RUO2, Rh02, OSO2 and 
Ir02 using first-principles calculations based on density functional theory (DFT). Firstly, 
their structural parameters are optimized by first-principles total energy calculations. Sec- 
ondly, the elastic constants are given by finite strain technique under the framework of 
linear response theory. Then the incompressibility of these TMO2S are discussed. Finally, 



the hardness and Debye temperature are estimated. 

II. COMPUTATIONAL DETAILS 

All calculations are performed by the CASTEP code using ab initio pseudopotentials 
based on DFT. For the exchange-correlation functional the generalized gradient approxima- 
tion (GGA), as given by Perdew, Becke, and Ernzerhof (PBE), is applied. All the possible 
structures are optimized by the BFGS algorithm which provides a fast way of finding the 
lowest energy structure and supports cell optimization. In the calculation, the interaction 
between the ions and the electrons is described by using Vanderbilt's supersoft pseudopoten- 
tial with the cutoff energy of 380 eV. The Monkhorst-Pack fc-point grid with a fine quality 
Brillouin zone sampling of 27rx0.04 A~^ is used in the calculations. In the geometrical 
optimization, all forces on atoms are converged to less than 0.002 eV/A, all the stress com- 
ponents are less than 0.02 GPa, and the tolerance in self-consistent field (SCF) calculation 
is 5.0x10^^ eV/atom. The relaxation of the internal degrees of freedom is allowed at each 
unit cell compression or expansion. The elastic constants of TMO2 are obtained by using the 
finite strain technique. From the full elastic constant tensor we determine the bulk modulus 
B, the shear modulus G, the Young's modulus, E, and Possion's ratio z/ according to the 
Voigt-Reuss-Hill (VRH) approximation. [lOt] 

III. RESULTS AND DISCUSSIONS 

The tetragonal rutile structure belongs to space group PA2/mnm. There are six atoms 
per unit cell. Two metal atoms locate at (0, 0, 0) sites and four oxygen atoms position at 
(u, u, 0) sites, where u is the internal parameter. The metal atoms are surrounded by six 
oxygen atoms at the corners of a slightly distorted octahedron, while the three metal atoms 
coordinating each of the oxygen atoms lie in a plane at the corners of a nearly equilateral 
triangle. Equilibrium volumes, crystal constants, and internal parameters are shown in Table 
I. It can be seen that all the structural parameters optimized here agree well with previous 
results. We find that there are very slight differences in volume and internal parameter [u) 
between TMO2S studied. 

The calculated elastic constants c^ are given in Table II. It is easy to see that these 



constants Cij satisfy the Born-Huang stability criteria, 

Cii > 0{i = l,3,4,6),cii > Ci2, 
cii + C33 > 2ci3, 2cii + C33 + 2ci2 + 4ci3 > 0, (1) 

suggesting that they are mechanically stable. From the elastic constants calculated above, 
we obtain the bulk modulus B, shear modulus G, Young's modulus E, and Possion's ratio 
u and list them in Table II, which are important to understand the elastic properties of 
TMO2S. For RUO2, our results agree well with the available theoretical and experimental 
values, indicating that PBE-GGA can be used to calculate the elastic properties of TMO2. 
Os02 and Ir02 have the close bulk moduli and Young's moduli, which are higher than those 
of Ru02 and Rh02. Higher bulk modulus indicates strong incompressibility. From Table II, 
we conclude that Ir02 and OSO2 are more difficult to be compressed than the two others, 
meaning that Ir02 and OSO2 are hard solids. The incompressibility of these TMO2 can 
also be seen from relative volume as a function of pressure (Fig. 1). While there isn't a 
large difference in the shear modulus at four TMO2S, it's roughly the same. Also we note 
that there are few difference (13 GPa at most) between c^A and shear modulus G for four 
dioxides. 

Employing the correlation between the shear modulus and Vickers hardness reported 



by Teter [12|, the hardness of four TMO2S is estimated. By using the calculated G values 
reported in Table II and the correlation given by Teter we find a theoretical Vickers hardness 
of approximately 14.1 GPa for RUO2 , 12.3 GPa for Rh02, 14.8 GPa for OsOa and 14.6 GPa 
for Ir02. It is well known that the hardness of superhard materials should be higher than 
40 GPa, so that TMO2 are not superhard materials. 

Possion's ratio refiects the stability of a crystal against shear. This ratio can formally 
take values between -1 and 0.5, which corresponds respectively to the lower limit where 
the material does not change its shape, and to the upper limit when the volume remains 
unchanged. All the calculated Poisson's ratios (Table II) are bigger than 0.25, which means 
that there are strong elastic anisotropy in these TMO2. In order to predict the brittle and 
ductile behavior of solids, Puch introduced the ratio of the bulk modulus to shear modulus 
of polycrystalline phases. A high (low) B/ G value is associated with ductility (brittleness) . 
The critical value which separates ductile and brittle materials is about 1.75. The bigger 
difference obtained between bulk and shear modulus indicates higher ductility of these TMO2 



(Table II). 

Now, we discuss the elastic anisotropy of TMO2. For the rutile structure, the directional 
bulk modulus along crystallographic axis and the percent elastic anisotropy defined by Chung 
and Buessem are discussed (Table III). From Table III, it can be seen that the directional 
bulk modulus is larger along the c axis than that along the a{b) axis, indicating that the 
incompressibility along the c axis is stronger than that along the a (b) axis. This anisotropy 
of compression along different lattice axis can be clearly seen from the increasing fractional 
axis compression c/a versus pressure (Fig. 2). On the other hand, the percentage of elastic 
anisotropy involves both the percentage anisotropy in compressibility Ab and in shear Ac- 
A value of zero represents elastic isotropy and a value of 1 (100%) is the largest possible 
anisotropy. The calculated Ab and Ac are listed in Table III. Obviously, these TMO2S 
possess strong shear anisotropy and weak bulk anisotropy. 

Further, the Debye temperature G^ is also discussed because the Debye temperature 

relates to many physical properties of materials, such as specific heat, dynamical properties, 

and melting temperature. It can be calculated from the average wave velocity Vm by the 
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where h is Plank's constant, k is Boltzmann's constant, N^ is Avogadro's number, p is den- 
sity, M is molecular weight and n is the number of atom in the molecule. Vm is approximately 
given by 
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where vi and Vt are the transverse and longitudinal elastic wave velocity of the polycrystalline 
material and are given by Navier's equation |l4| 
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The calculated density, longitudinal, transverse, and average elastic velocities and Debye 
temperatures Od for four TMO2S are given in Table IV. Debye temperatures 0d of RUO2 is 
comparative to the value of 866.4 K in ReB2 15| and ones of OSO2 and Ir02 is comparative 
to the value of 691 i^ in OsNa. 



IV. CONCLUSION 

In summary, the elastic properties of transition metal dioxides with rutile structure are 
investigated using first-principles calculations under the framework of density functional 
theory within generalized gradient density approximation. Elastic constants, bulk modulus, 
shear modulus, and Possion's ratio are given. Os02 and Ir02 have stronger incompressibility 
than Ru02 and Rh02. The estimated hardness is 14.1 GPa for Ru02, 12.3 GPa for Rh02, 
14.8 GPa for Os02 and 14.6 GPa for Ir02, which show that these dioxides are not superhard 
solids. Debye temperatures calculated by elastic constants and density are 846 K, 795 K, 
674 K, 674 K for RuOa, RhOa, OsOa, and IrOs, respectively. 
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FIGURE CAPTIONS 



Fig.l: V/Vq as a function of pressure for seven TMO2. 



Fig. 2: c/a ratio as a function of pressure for seven TMO2. 



TABLE CAPTIONS 

TABLE I: Equilibrium lattice parameters, Fo(A^), a (A), c (A), density p (g/cm^), and 
internal parameter u along with the available experimental values. 

TABLE II: Zero-pressure elastic constants Cij (GPa), the isotropic bulk modulus B (GPa), 
shear modulus G (GPa), Young's modulus E (GPa) and Possion's ratio u. 

TABLE III: The bulk modulus along the crystallographic axes a, and c {Ba and Be) for 
TMO2. Percent elastic anisotropy for shear and bulk moduli Aq (in %), Ab (in %) and 
compressibility anisotropy factors Ab^ and Ab^- Here, Ab^ =^=1, Ab^ —^■ 

TABLE IV: The longitudinal, transverse, average elastic wave velocity {vi, Vt, and Vm in 
m/s), and Debye temperature Qd {K) at the theoretical equilibrium volume. 



TABLE I: Li et al. 



Vq a c u Reference 

Ru02 63.1488 4.5076 3.1080 0.3055 This 

4.4919 3.1066 0.306 [1] 

Rh02 66.7652 4.6051 3.1482 0.3067 This 

Os02 64.3639 4.4943 3.1865 0.3078 This 

4.5003 3.1839 0.308 [1] 

Ir02 67.4532 4.5899 3.2018 0.3082 This 

4.4983 3.1544 0.307 [1] 

65.79 4.541 3.191 0.3085 [5] 



TABLE II: Li et al. 



cii C33 C44 Cm C12 cvi B G E v Ref. 
RUO2 303 536 117 221 240 190 260 104 275 0.3239 This 

299 558 114 227 246 199 258 [9] 

270 [llj 

Rh02 268 498 103 190 216 166 232 91 242 0.3263 This 
OSO2 337 631 109 243 270 210 292 109 291 0.3338 This 
Ir02 344 605 117 219 274 223 299 108 290 0.3381 This 



TABLE III: Li et al. 



Ba Be Ab, Ab Ag 
Ru02 633 1339 2.1155 1.66 21.01 
Rh02 560 1219 2.1751 1.89 22.24 
Os02 700 1573 2.2458 2.02 21.93 
Ir02 718 1600 2.2284 1.76 18.83 



TABLE IV: Li et al. 



P Vi Vt Vm Qd 

Ru02 6.9995 7548 3852 4316 846 
Rh02 6.7117 7261 3687 4133 795 
Os02 11.4687 6180 3086 3463 674 
Ir02 11.0413 6334 3133 3516 674 
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FIG. 1: Li et al. 
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FIG. 2: Li et al. 



